Transition from static to kinetic friction: Insights from a 2D model 
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We describe a 2D spring-block model for the transition from static to kinetic friction at an elastic 
slider/rigid substrate interface obeying a minimalistic friction law (Amontons-Coulomb). By using 
realistic boundary conditions, a number of previously unexplained experimental results on precursory 
micro-slip fronts are successfully reproduced. From the analysis of the interfacial stresses, we derive 
a prediction for the evolution of the precursor length as a function of the applied loads, as well as an 
approximate relationship between microscopic and macroscopic friction coefficients. We show that 
the stress build-up due to both elastic loading and micro-slip-related relaxations depend only weakly 
on the underlying shear crack propagation dynamics. Conversely, crack speed depends strongly on 
both the instantaneous stresses and the friction coefficients, through a non-trivial scaling parameter. 



Frictional interfaces are important in many areas of 
science and technology, including seismology [2j , biology 
[H, 0| and nanomechanics [E| . Whereas a satisfactory pic- 
ture of the steady sliding regime of such interfaces has 
been developed during the last twenty years [fHU, the 
dynamics of the transition from static to kinetic friction 
remains elusive. During the last decade, a renewed in- 
terest has grown in such transitions, due to experimen- 
tal studies that directly measured the local dynamics of 
frictional interfaces [9l-ll3j . They have shown that macro- 
scopic sliding occurs only after shear crack-like micro-slip 
fronts have spanned the entire contact interface. 

Experimentally, micro-slip front nucleation, propaga- 
tion and arrest was shown to be controlled by the in- 
stantaneous stress field at the interface. Fronts nucle- 
ate preferentially at the trailing edge of the contact area 
an effect explained cither by the enhanced 



shear stress near the loading point in side-driven systems 
f]~sl | or by a friction-induced pressure asym- 
metry in top-driven systems 0, [l9[ • Fronts can arise well 
below the macroscopic static friction threshold and arrest 



before the whole contact area has ruptured [14j - |16l ] . The 
length and number of these precursors depends on the 
precise way in which shear [14| and normal [16j forces are 
applied. Moreover, precursors are associated with signifi- 
cant changes in the spatial distribution of the real contact 
area [3], a quantity related to the local interfacial pres- 
sure. Finally, the propagation speed of micro-slip fronts, 
which covers a wide range [§- HI 2(|, correlates with the 
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local shear to normal stress ratio at nucleation 

Theoretically, some aspects of these observations have 
been studied using one-dimensional (ID) models. The 
conditions leading to a large range of front velocities were 
addressed using a ID spring-block model with a time- 
dependent friction law [la ]. The role of an asymmetric 
normal loading on the length of precursors was considered 
using a ID spring-block model with Amontons-Coulomb 
(A-C) friction and different normal forces ascribed to dif- 
ferent blocks [l6[ . The dependence of the series of pre- 
cursors on the friction-induced pressure asymmetry was 



described, for A-C friction, using a quasi-static ID model 
[l9| . In these three studies, the normal stress distribution 
was assumed to be uniform (lil ] or linear (either fixed (l6| 
or friction-dependent [3]). Such assumptions impede 
quantitative comparison with experiments since determi- 
nation of the actual stress field requires solution of the 
elastic problem for the two bodies in contact, including 
their precise geometry, elastic properties and boundary 
conditions, not only at the frictional interface but also 
on all their other boundaries. As a first step towards 
such complete description, wc present a minimal 1+1D 
(along + orthogonal to interface) model for the transi- 
tion from static to kinetic friction of an elastic slider on 
a rigid substrate. We show that this model, by enabling 
realistic boundary conditions, is sufficient to reproduce a 
series of still unexplained experimental observations. 

We used the 2D spring-block model sketched in Fig. 
[Ha). The slider has mass M and sizes L and H in the 
horizontal (x) and vertical (z) directions, respectively. It 
is divided along a square lattice into N = N X N Z blocks 
of mass m = M/N. Blocks are coupled to their four 
nearest neighbours and their four ncxt-ncarcst neigh- 
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FIG. 1. (a): Sketch of the 2D spring-block model, (b) Blue: 
Typical loading curve, from static contact to macroscopic 
stick-slip, for setup 1 (/i=5mm, L=140mm, _Fjv=3500N). 
Green: Locations of slipping regions for all micro-slip event. 
(I) corresponds to the trailing (leading) edge. 
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bours by springs of equilibrium lengths I — L/(N X — 1) = 
H/(N Z — 1) and \/2l and stiffnesses k and fc/2, respec- 
tively, giving an isotropic elastic model with Poisson's 
ratio f/3. The force exerted on block i by block j is 



thus kijfrij 



I ■■)'- 



when blocks are connected, oth- 
erwise, where x = (x,z), Ax,j = Xj — Xj, ry = |Ax,;j| 
and kij and Zy are the stiffness and equilibrium length of 
the spring connecting blocks i and j. Block oscillations 
are damped using a viscous force rj(xj — x,-) on the rela- 
tive motion of connected blocks. We chose the coefficient 
i] = vO-1 km so that blocks are underdamped and event- 
triggered oscillations die out well before the next event. 
All r\ satisfying these conditions gave similar results. 

Boundary conditions were designed to model two dif- 
ferent experiments described in the literature: [3 (setup 
f) and (setup 2). They differ in the way the top 
blocks are loaded. For setup f , they are glued to a rigid 
rod of mass 75. 6g, itself submitted to a normal force 
Fjy and coupled to a soft "spring mattress" of stiffness 
0.4MN/m, the effect of which is modeled with a restoring 
torque proportional to the rod's tilt angle. For setup 2, 
the top blocks are submitted to a linear time-independent 
distribution of vertical forces jf- (l + 2t jv jv ;~ 1 6>^) , where 

6 € [—1,1] controls the pressure asymmetry. In both 
setups, the bottom blocks lie on an elastic foundation 
of modulus kf = k, i.e. each block is submitted to a 
vertical force of amplitude pi = kf \zi\ if z.; < or oth- 
erwise, where Zi is the vertical displacement of block i. 
All kf > k gave similar results, so that the substrate can 
be considered rigid compared to the slider. Both verti- 
cal boundaries are free, except for a horizontal driving 
force Ft — K(Vt — Xh) applied on the left-side block 
situated at height h above the interface, where Xh is the 
x-displacement of this block. This models a pushing de- 
vice of stiffness K driven at a small constant velocity V. 
The amplitudes fi of the friction forces on the bottom 
blocks follow the minimalistic local A-C friction law with 
static and kinetic friction coefficients fi s and [ik < fJ-s- 
If ii = 0, fi balances all horizontal forces on block i up 
to n s pi, at which slip initiates; then fi = ±[ikPi when 
ii 0. The 2A^ equations of motion are solved simulta- 
neously using a fourth order Runge-Kutta integrator on 
a uniform temporal grid of resolution At. 

Using this model, we simulated the transition from 
static to kinetic friction for the various loading config- 



urations reported in 14J and [161 ] . Model parameters, 
chosen in accordance with experiments (2l| , are given in 
the caption of Fig. [2] Figure \V[b) (blue) shows a typi- 
cal loading curve for setup 1, from static contact up to 
macroscopic stick-slip. As in ID models, and in agree- 
ment with experimental results, we observe, well before 
macroscopic sliding, a series of partial force relaxation 
events, corresponding to precursors, which all nucleate 
at the trailing edge, and extend over increasing lengths 
L P (Fig. [Hb), green). 
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FIG. 2. Normalized precursor length L p /L as a function of 
normalized shear force Ft/Fn- (a) Setup 1 for various config- 
urations: L=140 or 200mm, h ranges between 2.5 and 15mm, 
Fjv=1750, 2700 or 3500N, #=75mm, _K=4MN/m, M=75.6g, 
V=0.7875mm/s, Z=2.5mm, fc=13.5MN/m, ^ s =0.7, Mfc=0.45, 
At=0.2[is. (b) Setup 2 for 0=0.833 (+), (□) or - 
0.833 (V), L=100mm, #=20mm, if =0.8MN/m, M=12g, 
V=0.45mm/s, £V=400N, h=2mm, !=lmm, fc=9.375MN/m, 
/x s =0.7, fj,k=0A5, At=0.09u.s. Vertical dashes: [it- Solid 
lines: prediction obtained with [i 3 = fit- 



We first focus on the dependence of L p on the applied 
tangential force Ft just after relaxation. We simulated, 
for setup 1, different slider lengths L, pushing heights 
h and normal forces F]y. The behaviour under different 
conditions differ by the number of precursors occurring 
along the transition: the increase in both L p and Ft be- 
tween precursors scales almost linearly with h for small 
h, so that larger L and smaller h yield more precursors. 
However, Fig. [U^a) shows that the results for all condi- 
tions can be collapsed on a single curve by plotting L p / L 
as a function of Ft/Fn- The same collapse was found in 
[l4| , with a very similar shape for the non- linear increase 
of L p /L with Ft/Fn. In particular, we reproduce the 
transition from a roughly linear increase up to L p /L ~0.5 
to a more rapid growth for longer precursors. We empha- 
size that a ID model with homogeneous normal loading 
would produce a purely linear increase. In Fig. |^b) we 
compare, for setup 2, the evolutions of L p / L as a function 
of Ft I Fn for three different linearly asymmetric normal 
loadings of the slider. Qualitatively, the lower the normal 
load on the trailing edge, the lower the threshold force re- 
quired to nucleate precursors, and therefore the lower the 
tangential force Ft necessary for the precursor to reach 
a given length L p , which explains the relative positions 
of the three curves in Fig. [UJb). Again, the simulated 
curves are in excellent agreement with the experimental 
results in Fig. 8 in la ], much better than the ID sim- 
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ulation (Fig. 

both setups clearly demonstrate that, by enabling realis- 
tic boundary conditions, 2D models do offer significantly 
improved agreement with experiments. 

We now use our simulation results to gain insight into 
how shear cracks nucleate, propagate and arrest at an 
interface obeying A-C friction. In particular, we will ex- 
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FIG. 3. Spatial distribution of stresses at three stages of 
the transition. Dotted lines: initial configuration (^=0). 
Dashed (solid) lines: at crack initiation (arrest) for the 9th 
precursor shown in Fig. [ljb). Blue: r. Black: (x s p. Red: fikP- 
Square (disk): nucleation (arrest) point, (b) Spatial distribu- 
tion of p, normalized by the initial distribution, after each of 
the precursors shown in Fig. [lib). Inset: Similar distribution 
for system-sized events during macroscopic stick-slip. 



amine the role played by the evolution of the shear and 
normal stress distributions at the interface. We empha- 
size that, in 2D, the latter is a result of the simulation, 
not an assumption, which is required in ID. In Fig. G3a) 
we show the typical stress evolution for setup 1. Ini- 
tially (Ft = 0) the normal stress p(x) is symmetric, with 
edge effects related to the flat punch geometry of the 
contact. The shear stress t(x) is antisymmetric, due to 
friction-frustrated Poisson expansion. These stresses are 
in excellent agreement with those expected from contact 
mechanics and those measured in [17(. Application of a 
tangential force at the trailing edge modifies slightly the 
normal stress field and modifies significantly the shear 
stress field, with a large increase near the trailing edge 
over a distance of order h. The local slipping threshold 
t(x) = fi s p(x) is therefore reached first near the trail- 
ing edge and a micro-slip front corresponding to the first 
precursor nucleates there. It then stops after propaga- 
tion over a finite distance L p and the whole scenario is 
repeated until the leading edge is reached. Figure 03a) 
illustrates this scenario for a typical precursor event. 

We find that the shear stress r(x) just after a precur- 
sor is always very close to /ikp(x) over the whole slipped 
length x < L p . This shows that the arrest state of the 
interface is only weakly dependent of the static friction 
coefficient fi s , and is controlled primarily by the kinetic 
friction coefficient fik ■ We emphasize that this behaviour 
is not specific to A-C friction, but remained true for slip- 
weakening friction, provided the weakening distance is 
smaller than a few /zm. Based on this robust behaviour, 
we propose the following procedure to predict, for given 
loading conditions and a given the non-linear evolu- 
tion of L p with Ft'- We run the model for equal friction 
coefficients (/x s = /i/.). The transition to kinetic friction 
in this simplified model is smooth, with a continuously 
growing micro-slip region the length of which is shown 



as a function of the applied force in Fig. [2] (solid lines). 
For both setups and for all loading configurations, this 
curve is in good quantitative agreement with the curve for 
the length of precursors as a function of the arrest force. 
Such agreement is due to the fact that any arrest state in 
the full model is very similar to the state reached in the 
simplified model for the same force Ft', (i) in the micro- 
slip region t(x) = Hkp(x) and (ii) in the stuck region 
t(x) arises from almost identical boundary conditions. 
Because internal dynamics are fast compared to changes 
in external loading, the prediction is ^-independent and 
can be obtained using simple equilibrium calculations, 
i.e. it does not require a complete dynamical simulation. 

The very last precursor in the full model propagated 
over almost the entire interface and left a shear stress 
that was equal almost everywhere to HkP, yielding a total 
shear force Ft — /Ufc-Fiv- The last increment of tangen- 
tial force required to trigger the first system-sized event 
brings the vicinity of the trailing edge to its threshold, 
whereas the shear stress on the rest of the interface is 
essentially unchanged, yielding a maximum total shear 
force only slightly above fik^N (see Fig. QJb)). This 
maximum force is classically interpreted as /x™ acro FAr 
with ^™ acro the macroscopic static friction coefficient. 
Therefore our results suggest that, in side-driven systems, 
whatever the value of p, s , ^™ acro ~ The difference 
between the macro- and microscopic static friction coef- 
ficients, already discussed in recent ID models flU [l9j. 
provides a possible explanation for the anomalously high 
values of fj, a compared to ^™ acro reported in [T^] ■ 

Not only the shear stress r but also the normal stress 
p is altered along the transition. Figure EJb) shows p(x), 
normalized by the initial distribution (at Ft = 0), af- 
ter each of the successive precursors. The normal stress 
is found to be significantly decreased along the path of 
the precursor that just occurred, by up to 30% whereas, 
apart from edge effects, it is mainly unchanged in front 
of it. The normal stress then assumes a reproducible dis- 
tribution in the macroscopic stick-slip regime (Fig. EJb), 
inset). Recalling that, at normally loaded rough contact 
interfaces, normal stress is robustly found to be locally 
proportional to the area of real contact, the results of Fig. 
EJb) show very good agreement with measurements of the 
real area of contact reported in [3] (Fig. 4a of Di- 
rect quantitative comparison is not possible, mainly be- 
cause normal stress has a constant integral (normal force 
Fn imposed) whereas the total real area of contact is 
not a conserved quantity, but typically decreases by 20% 
across the transition 0, (HJ . 

Let us now focus on the propagation dynamics of the 
simulated frictional shear cracks. Figure Efla) shows the 
evolution of the local crack speed v c as a function of posi- 
tion x along the interface for a typical precursor event in 
setup 1. v c (x) decreases from around 2400m/s near the 
trailing edge to around 400m/s just before crack arrest 
and appears to be strongly correlated to the shear to nor- 
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FIG. 4. (a) Local crack speed v c (solid line, spatial average on 
5 neighbouring blocks) and local shear to normal stress ratio 
r/p (dashed line) as a function of crack length x for the 11th 
precursor event shown in Fig. QJb). (b) Inset: v c as a function 
of r/p. Runs: 1 (left), 2 (center), 3 and 4 (right). Main: v c as 
a function of ll£=Oh . R uns: \ ( V ), 2 (o), 3 (A), 4 (■)• The 
longitudinal (shear) wave speed is 1677m/s (968m/s). 
Run 1: L=140mm, Fjv=3500N, /i=5mm, ^ s =0.5, jij—0.2. 
Run 2: L=140mm, fW=3500N, /i=5mm, /i s =0.55, ^fc=0.3. 
Run 3: L=140mm, Fjv=3500N, ft=7.5mm, /i s =0.7, /x fc =0.45. 
Run 4: L=200mm, F]v=2700N, /i=10mm, /x s =0.7, ^=0.45. 



mal stress ratio r(x)/p(x) (Fig. HJa)). In fact, when v c 
is plotted as a function of r/p, all points corresponding 
to all locations x along all successive events in a given 
simulation collapse on a single curve (Fig. Sl^b), inset). 
Different geometries or loading conditions yield the very 
same curve but different friction coefficients yield differ- 
ent curves (Fig. 0Jb), inset). All these curves can then be 
collapsed on the same master curve when v c (x) is plotted 
as a function of the non-trivial parameter r ^|/ P l'^ fc ~ Mfc , 
which represents the local distance to the slipping thresh- 
old and extends existing parameters @, [23| to spatial 
heterogenities in normal stress p(x). These results are 
a generalization, accounting for any value of the micro- 
scopic friction coefficients, of a similar collapse obtained 
for experimental data, from system-sized events only, us- 
ing r/p as a parameter (Fig. 3 in [l3l)- Our results also 
suggest that, if probed, precursors would follow the same 
experimental curve as system-sized events. 

In surprising contrast with the excellent agreement 
found with experiments up to now, the master curve of 
Fig. BJb) exhibits strong discrepancies with its experi- 
mental counterpart. The shape as well as the explored 
ranges of both v c and r/p are different. In particular, 
our model does not produce very slow micro-slip fronts 
like those observed in various experiments 17J, [20| . Most 
likely this is because the minimalistic friction used here 
(A-C) lacks some time-dependent ingredient necessary to 
yield slow fronts, like those in [l8j |. 

All our results suggest two distinct levels of description 
of the transition from static to kinetic friction. First, a 
kinematic description of (i) the slow evolution of inter- 
facial stresses between events and (ii) the stress condi- 



tions at crack nuclcation and arrest. Second, a dynamic 
description of the fast propagation of micro-slip fronts 
along the interface. Crack dynamics appear to depend 
crucially, via the friction law, on the kinematic stresses 
at crack initiation. Conversely, the kinematic descrip- 
tion was found essentially independent of the underlying 
dynamics: we could successfully reproduce all available 
experimental kinematic results, even with an unrealistic 
friction law. The key observation yielding this surpris- 
ing success is that, at crack arrest, t(x) = fj,i-p(x) over 
the whole slipped region. Since the other 2D boundary 
conditions on the slider could be accurately taken into 
account, our results strongly suggest that a very simi- 
lar arrest condition holds in the experiments. We thus 
believe that any friction law leading to such arrest state 
will produce as good kinematic agreement as A-C fric- 
tion, irrespective of the dynamic way of reaching this 
arrest state. 
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